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The bilinear-biquadratic spin 1 quantum chain 

L 

H = J2 coses n ■ S n+ i + sm0(S„ • S„+i) 2 (1) 

n=l 

has been the object of very intense investigation in the last twenty years. One of the goals has been 
to determine the properties of the ground state and the nature of the low-lying excitations in the 
thermodynamic limit. In particular, whether they form a continuum with the ground state (critical 
chain) or whether they are separated by a finite gap (massive chain). Haldane's prediction that 
the Heisenberg antiferromagnct (9 = 0) should be massive for integer spin has been thouroughly 
verified by a variety of numerical methods The AKLT chain (tan# =5) has been rigorously 
proven to have a valence bond ground state and a nonzero gap Q . Other cases of (||) have been 
analyzed in depth due to their integrability. They are the critical SU (3)- invariant Sutherland spin 
chain (8 = 7r/4), solved by nested Bethe-ansatz Q; the Babujian-Takhtajian chain (9 = — 7r/4), 
also solvable by Bethe-ansatz and also gapless |||| . The third case is the purely biquadratic spin 
chain (9 = —ir/2), which, after some controversy, has been found to have a gap A ~ 0.1731788 [Q; 
an identical gap has been revealed in the same chain, but this time with free boundary conditions 

L-l 

HbQ = - yi(Sn • S„ + i) 2 (2) 

71=1 

and L even ||. 

Now, at least for free boundary conditions, it can be shown that (|J) has a band of gapless 
excitations when L is odd and, of course, L — > 00. This points to the fact that antiferromagnetic 
quantum spin chains can in principle have very different properties according to the parity of their 
length. 

To prove this statement, one observes that the spectrum of (|j) can be mapped ^ into the 
spectrum of the s/(2) g -invariant open antiferromagnetic XXZ chain with a boundary field 

1 L_1 sinh 
Hxxz = J2«<+l + " cosh 7 « +1 ) + — 2-H<r? - of) (3) 

71 = 1 

The equivalence holds when cosh 7 = | , sinh 7 = At this point both (Q) and (J3j> can be written 
as sums over the generators {e n }^z\ of the Temperley-Lieb algebra ||,|| 

e 2 = 3e„ e„e„±ie„ = e„ [e„, e n >] = \n-n'\>2 (4) 

H bQ = H X xz- \{L-l) (5) 

Eq. (|5|) must be understood as a statement about the spectra. The two representations of ([!]), 
appearing in (Q) and (||), are of course different since the two Hamiltonians live in Hilbert spaces 
of different dimensions, but their spectra should differ only in their multiplicities p|JTo|] (further 
comments on this point will be given at the end). Similar mappings, via Temperley-Lieb algebra, 
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have shown to be very fruitful in studying statistical mechanics models and related quantum spin 



chains, one notable example being the Potts model pO 11|. 

The spin chain (||) has been solved by means of the coordinate Bethe-ansatz (BA) [Q. In fact, 
it was relying on the Temperley-Lieb mapping and the numerical as well as analytical solution of 
the relevant BA equations that the ground state energy and the gap of (||) were computed in 
As in all BA solvable systems, the spectrum of (|^) is expressed in terms of rapidities {eti, . . . a n } 
that solve a set of coupled equations. In the case at hand, they are 

11™ I 

E (e(oi-a*;7) + e(aj+Ofc57)) = f j = l...n (6) 

sinh(a; + f ) 



0(a; x) d = —i In 



sinh(a; - f ) 



(X 

2 arctan(tan — coth x) (7) 



Solutions are labelled by the set of positive integers {ij}" = i. The branch cut in (^) is chosen 
to make Q(a;x) a differentiable, increasing function for a real, with 0(— a;x) = —Q(a,x) and 
Q(a + 2ir; x) — Q(a; x) + 2ir. Only positive real rapidities in the (0, it) interval need to be considered 
here. Given a solution of (||), the corresponding energy and spin Sf ot — \ J2j=i a j arc @ 

1 " L 

E=-{L-l) cosh 7 - 2sinh 7 ^ e'(a,; j) S* tot = -^-n 

i=i 

It's now crucial to determine the range of allowed vacancies for the set {Jj}™ =1 . This is done by 
rewriting (^|) in terms of the counting function |D3] , defined as 



de f 7n , 1 



7T v ' 2 ' 2nL 



Z L (a) = J -6(a;^) + — 6(a; 7 ) + 6(2a; 7 ) - ^ > J e(a-a fc ; 7 ) 



1 



2ttL ^ 

fe=-T 



z L{ctj) = j- j = -n, -n+ 1, .. .,n- l,n (8) 

Rapidities have been doubled by reflection through a — fl^] : {ai,...,a„ ; > 0} — > 
. . . , a_i, 0, ai, . . . , a„ | a_j = — a^}, so that (||) is completely equivalent to (^|). Since 
the rapidity range is (0,7r), or (— 7r,7r) after reflection, and the counting function is monotonically 
increasing, the largest allowed Ij is read from LZl(t:), or, conceptually more correctly but with 
equal numerical result 

lim LZ L {a 3 ) =L-n+l (9) 

aj — >7r 

The limit does not depend on the numerical value of the remaining at, k ^ j. Eq. (||) would 
suggest I max — L — n + 1. Actually this value is forbidden because the wave-function vanishes 
when one rapidity is ir (a is related to the variable k of (l2) by e lk = (e la — e~ 7 )/(l — e la ')). 
Hence the largest possible integer is 

Imax = L-n (10) 
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If L is even, n — £ for the ground state (sector S% ot = 0) and there's exactly h vacancies for h 
positive integers. But if L is odd, the lowest energy states have n — (sector Sf ot = 1/2) and 
there are vacancies for integers, hence one hole. Altogether a band of configurations, 
specified by a sequence of closely packed integers with one hole JW, 1 < < ^±1. The ground 
state for finite L has 1^ = -^yl , that is at the edge of the band, as it will be proven shortly. 

Denote by {<x,}" =1 the rapidities of the state with jw = by {a^}™ =1 the rapidities of 

any other state in the band (1 < 1^ < ^-^-) and by zf\a) and Z^\a) the relevant counting 
functions. One defines a hole rapidity a^) by 



4V°) = V 



T 



This number does not belong to the set {a^}" =1 , but it is convenient to include it, adding and 
subtracting its contribution to the counting function, which reads after reflection 

n+1 



4%) 



1 T 1 



e(a; 7 ) + e(2a; 7 ) 



1 

2nL 



J2 ©(a -4; 7) 



k——n— 1 



2~T ( e(a ~ a(M ; 7) + 0(a + aW ; 7) ) 



so that 



Z^(a k )-Z^(a' k )=Q 



L — l , L 
— - — < k < — 



(11) 



For these values of k, define 5a k = a' k — a k = 0(1/ L). The two (after reflection) unpaired edge 
rapidities are handled separately. Since {oj }™ = _ n fill the interval (— it, w) in the thermodynamic 



limit m& . one can set a'. 



o(l) and a 1 



o(l) as L — > oo. The next steps follow 



a well-established method for handling the BA equations |lj,|l5|. Defining Act* — cxk+i — ct k 



0(1/ L) and the backflow 8(a) = lim^^ 
cancel and the terms 0(1 /L) yield 



0(1), it is found that the terms 0(1) in (jl 



5(a) + — 
Zir 



d/3G'(a - /?; >y)5(0) = — Q(a- vr; 7 ) + 6(q + vr; 7 ) 
Zir L 



-Q(a - a (h) ; 7) - 0(a + a {h) - 7) 

A£(c>)) = lim (^(qW) = -sinh 7 / daQ' 

+2 S inh 7 (e'(aW; 2) - 0'^ ~)) 



(a; -)5(a) 



(12) 



(13) 



It is straightforward to solve (|l2] ) by Fourier transform. It is actually better to reduce it to an 
equation for 5' (a) by differentiating and then using 5(ir) — S(—n) which follows from ( [l2] ) itself. 
5' (a) is substituted after (Ol) has been integrated by parts. The result is 



AE(a^) =e(aW)-e(7r) 



< a {h) < 7T 



(14) 
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e( a ) = 2 S inh 7 g ^^ = 2sinh7 ffl dn( ^ ;fc) 
^— 4 2 cosh 771 tt tt 

71— — OO 

K'(k) _ 7 
K(k) ~ 7T 

where K(K') is the real (imaginary) quarter period and k the modulus of the elliptic dn function 
fL7f . This shows that AE(a^) > and it goes to when — > 7r, as claimed. Moreover, the 
whole band found for L odd in S* ot = 1/2 is degenerate with a twin band in the sector S^ ot = —1/2. 
In fact, despite the boundary field, spectra in the sectors S% ot = S and S% ot = — S coincide because 
of the symmetry implemented by the unitary operator J[]^ =1 cr^, where Ia n I = <7L+i-n- 

Excitations of the XXZ chain are sometimes called spinons. When L is even they appear in 
even number |l3|JlS|] , hence they are assigned a spin S z = 1/2. Eq.(jl4|) shows that, when L is 
odd, the ground state sector must contain one spinon which can be in a whole band of dynamical 
states with varying energy. Such states cannot be labelled by a linear momentum because (|J) is 
not traslationally invariant. 

The case L even is known. In the sector S% ot = 1 there are n = j — 1 rapidities; this, from ([To|), 
implies 4 + 1 vacancies for ^ — 1 integers, hence two holes leading to 

A£( tt (W), a (M) ) = e(a (M) ) + e(a ( "' 2) ) < < o^< 2 > < it 

This coincides with the well know result for the periodic XXZ chain, the only difference being 
that now hole rapidities lie in the restricted range (0, tt). At cosh 7 = § the gap is of course 
2e(7r) ~ 0.1731788 §. 

Other quantities have a different limit for L odd or even, one example being the surface energy 
e( s \ defined by 

£ (0) (L) = e L + + o(l) L -> 00 

For L even it has been calculated in |l£j] through an analysis of finite size corrections. Since the 
XXZ periodic chain has the same eo = cos 2 h 7 — sinh7(l + ^Y^=i i+l^-i ) lHH but no surface 
energy, can also be found from 

2£ (0) (L) - E {0 ' p) (2L) = 2e (s) + o(l) L -> 00 

where the label "p" stands for periodic boundary conditions. This is an order 1 calculation which 
bypasses the intricacies of finite size calculations. For the periodic chain of length 2L, the ground 
state has n = L (HQ 

Z^\a) = -e(a; 2) - ^ £e( a - a k] 7) 

71 ~ 71 k=l 
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E^- p) (2L) = L cosh 7 - 2 sinh 7 ^ Q'(«fc! ~) 

fc=i 

Consider, when i is odd, the lowest-lying state in the previously discussed band of the open spin 
chain. Its (L — l)/2 rapidities, which become exactly L after reflection, will now be denoted 
{^•^i—fi^i) /21 whereas /2 wn ^ be ^ ne symmetrically distributed ground state 

rapidities of the periodic chain of length 2L. Hence 

17(C)) / / \ y(Oj>), \ j L-l £-l 

Z L ( a j) = Z 2L = 75 — < J < ~ ^~ 

The relevant equations in the thermodynamic limit are now 

5(a) + ^ / d/30'(a - /3; 7)^) = [©(«; 7) + ©(2a; 7) 



2e« 



/7T 
da9"(a; ~)<5(a) 
-7T ^ 



for the backflow (5(a) = lini[,_ >c>0 (as before Saj = a'j — aj, Aaj — oy+i — aj). The upshot is 
(sodd) cosh 7 . , / +^2-e- 2fc ^ 2 \ 



cosh2fc7 ^— ^ cosh(2fc + 1)7 



For L even, the calculation is slightly more involved, because now the numbers Ij are half-odd for 
the periodic chain of length 2L, yet it can be carried out along the same lines, leading to 

cosh 7 , +2?l- e - 2 ^ 



e (^-)=_^iI+2 S inh7X: 



cosh2fc7 
fc=i ' 

A similar method was used in |2lJ] for general boundary fields but, again, for L even only. When 
COSI17 = 3/2, the biquadratic chain surface energy is, from (j^), e^Q = e( s ) + 7/4, or 

e (sMd) _ 1J415986 e^ ven) ~ 1.6550092 

The second value coincides with that in . 

I have checked numerically, for small chains of odd length, that: 1. eq. (||) admit a solution for 
any value of within the allowed range; 2. the corresponding energies show up in the spectrum 
of (H|), as found by direct diagonalization, and = yields the ground state energy; 3. the 
same eigenvalues show up in the spectrum of also obtained by direct diagonalization. The last 
check was carried out in @] too (perhaps only for L even) . 

Actually, the eigenvalue problem of (||) has a direct BA solution that does not rely on the 
mapping (|o) , but the relevant BA equations do not seem to have been studied in detail. This 
would provide a further check. Finally, it should be pointed out that, very recently, peculiarities 
of L odd have been noticed also for the periodic XX Z chain in the critical regime 0]. 
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